Non-Ohmic critical fluctuation Hall conductivity of layered superconductors in strong 

electric fields 
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The excess Hall conductivity, resulting from thermal fluctuations of the superconducting order 
parameter, is calculated for a layered superconductor for an arbitrarily strong in-plane electric field 
and a perpendicular magnetic field in the frame of the time-dependent Ginzburg-Landau theory. The 
fluctuation Hall conductivity is suppressed in high electric fields much stronger than the longitudinal 
one. For high-temperature superconductors we predict a pronounced non-Ohmic behavior of the 
excess Hall effect near the critical temperature in moderate magnetic fields and electric fields of the 
order of 100 V/cm. 



The non-Ohmic behavior of the fluctuation conductiv- 
ity in strong electric fields, studied for the first time for 
the isotropic case in connection with the low temperature 
superconductors^*^ can be summarized by saying that 
reasonably high values of the electric field can accelerate 
the fluctuating paired carriers to the depairing current, 
and thus, suppress the lifetime of the fluctuations, which 
leads to a deviation from Ohm's law. For a layered su- 
perconductor, a situation very much resembling the crys- 
tal structure of the high temperature superconductors 
(HTSC) , the issue has been addressed theoretically start- 
ing from a microscopic approach'^ and subsequently in the 
frame of the time-dependent Ginzburg-Landau (TDGL) 
theory, in the Gaussian^ as well as in the self-consistent 
Hartree approximation.^ The non-Ohmic fluctuation con- 
ductivity was studied recently also in the presence of a 
perpendicular magnetic field^^ in the Hartree approxima- 
tion for the TDGL theory, revealing that the simultane- 
ous application of electric and magnetic fields results only 
in a slight additional suppression of the superconducting 
fiuctuations, compared to the case when the fields are 
applied individually. 

However, a treatment of a possible non-linearity of 
the off-diagonal (Hall) components of the fluctuation 
magnetoconductivity tensor in a high electric field has 
been lacking yet. The non- vanishing Hall current 
due to fluctuating Cooper pairs, ascribed to a hole- 
particle asymmetry^ and a complex relaxation time in 
the TDGL theory, was first calculated in the Gaus- 
sian approximation^ and subsequently improved by in- 
corporating the fluctuation interaction in the frame of 
a Hartree approach,* or based on the single particle- 
hole renormalization.^ All the above models for the so- 
called excess Hall conductivity represent however linear 
response approximations with respect to the longitudi- 
nal electric field and are therefore valid only for small 
magnitudes of the latter. 

Experimentally, several investigationsiSiiiiiSii^ of the 
Hall effect at high current densities up to 10^ Acm~^ have 
been reported, but with the aim of overcoming the vortex 
pinning and hence of testing its influence on the Hall 
conductivity. The intrinsic non-Ohmic effect was neither 
envisaged nor fortuitously evidenced, since it needs, as 
we shall see later, even higher current densities in order 



to be unambiguously discerned. 

In this paper we shall treat, in the self-consistent 
Hartree approach, the thermal fluctuation Hall conduc- 
tivity for a layered superconductor in a perpendicular 
magnetic field and for an arbitrarily strong in-plane elec- 
tric field. To our present knowledge this topic has not 
been treated yet, not even for the simpler cases of an 
isotropic superconductor, or in the Gaussian approxi- 
mation. For our purpose, we shall adopt the Langevin 
approach to the TDGL equation. ■'^-^ Keeping the same 
notations as in Ref. |^ the gauge-invariant relaxational 
TDGL equation governing the critical dynamics of the 
superconducting order parameter in the ^-th supercon- 
ducting plane writes: 
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Here m and rric are effective Cooper pair masses in the ab- 
plane and along the c-axis, respectively, s is the distance 
between superconducting planes, and e the elementary 
electric charge. The order parameter has the same phys- 
ical dimension as in the three-dimensional case and SI 
units are used. The perpendicular magnetic field B is 
generated by the vector potential A — (0, xB, 0), with x 
and y the in-plane coordinates, and the magnetization is 
neglected. The GL potential a = aoe is parameterized 
by ao = r?l2mil = ^'^/'^rn^Qc and e = In (T/Tq), with 
To the mean-field transition temperature, while ^0 and 
^Oc are the in-plane and out-of-plane coherence lengths 
extrapolated at T = 0, respectively. The real part of 
the relaxation time in the TDGL equation is giveni^ by 
Fq"^ = nh^/16m£^QkBT, while the imaginary part Fq ^77 
must be introduced in order to break the particle-hole 
symmetry and allow for a non- vanishing Hall current iSiS 
The Langevin white-noise forces Q (x, t) are correlated 
through (0 (x,t)Cp (x',<')) - 2r^^kBT6ix - x')S{t - 
t')^u'/s, where 6{x — x') is the 2-dimensional delta- 
function concerning the in-plane coordinates. The elec- 
tric field E is assumed along the a;-axis, generated by the 
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scalar potential = —Ex. In the chosen gauge, the fluc- 
tuation current density along the y direction in the Z-th 
plane, averaged with respect to the noise, writes: 

ieh , 
m 

m \ 

so that the fluctuation Hall conductivity is given by 

As mentioned, the quartic term in the thermo- 
dynamical potential will be treated in the Hartree 
approximationjSii^ which results in a linear problem with 
a modifled (renormalized) reduced temperature e = e + 

b(\^if)/ao. 

Following the same procedure as in Ref. 0, we intro- 
duce the Fourier transform with respect to the in-plane 
coordinate y, the layer index I, and time respectively, 
and also the Landau level (LL) representation with re- 
spect to the a;-dependence, through the relation: 

f dk dq f dto ^ , , , , 

M^,y,t)^J - j_ -j -^UnM (3) 
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where the functions u„ (a;) with n G N build the or- 
thonormal eigenfunction system of the harmonic oscil- 
lator hamiltonian, so that [—fPd^ + 4e'^B^x'^)un{x) = 
2heB{2n + l)u.n{x). 

Equation |^ can be evaluated into the matrix form, 
after applying the expansion 0: 



ao 



^ (M + 77P)„„, ^,in', fc, uj) = C,in, k, u;) (4) 



where the new noise terms C,q{n,k,Lo), corresponding 
to the expansion rule (0), are delta-correlated such as 
(Cgi (ni,fci,wi)C*, {n2,k2.L02)) = 2r^ifcBT(27r)35(fci - 
k2)5{qi — q2)6{ijJi — uj2)5nin2^ 8.nd where the two dimen- 
sionless symmetrical tridiagonal matrices M and P have 
the elements: 
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Here we have introduced the new variables: 
Ek 



p-i 
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and the reduced field magnitudes: 
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with Eq = IGVSfcsT / Tre^o defined as in Refs. iandU 

By solving Eq. Q), and taking into account the ex- 
pansion form (O, one obtains the correlation function of 
the order parameter: 



{ipi {x, y, t) i})* {x, y', t)) = —^—h 



(8) 
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■u,, {x - k') f (M + r/P)+ • (M + 77P)1 (g', J) , 



where the notation [•••]„„/ is to be understood as the el- 
ement of the inverted matrix. 

Before proceeding further we point out that the sums 
over the LL in Eqs. (O, Q and © must be cut-off at 
some index iVc, reflecting the inherent UV divergence of 
the Ginzburg-Landau theory. The classical^-'i^. procedure 
is to suppress the short wavelength fluctuating modes 
through a momentum (or, equivalently, kinetic energy) 
cut-off condition, which, in terms of the LL represen- 
tation writes*'^^ (heoB/m) {n + ^) < cap = ch?/2m£,'^, 
with the cut-off parameter c of the order of unity. A total 
energy cut-off was also recently proposed, whose phys- 
ical meaning was shown to follow from the uncertainty 
principle. However, in the critical fluctuation region the 
two cut-off conditions almost coincide quantitatively, due 
to the low reduced-temperature e with respect to c, so 
that we shall apply for simplicity the momentum cut-off 
procedure. In terms of the reduced magnetic field h, it 
writes thus h (iVc + 5) = c/2. In this way, the matrices 
M and P are truncated at iVc + 1 lines and columns. 

The value of r] can be inferred from the microscop- 
ical theory if one considers the energy derivative A/"' 
of the density of states Af at the Fermi level ep, and 
it writesS 77 = — (kBT/ep) ce, where the parameter a 
amounts in the BCS model to a = AepM' /ngBCS-^"^ , 
with gncs (> 0) the BCS coupling constantii^ Since 
ep is for HTSC of the order of 10^ K (in ke units) 
and the hole-particle asymmetry parameter a, inferred 
from fits of excess Hall effect data^^-^- with the mod- 
els from Refs. and 0, turns out to be of the order of 
10^2 10^^, we conclude that 77 is a small parameter, 
reflecting also the small Hall angle. We shall therefore 
expand the inverted matrix in Eq. ^ only up to the 



linear term in 77, such as 

Q 
Q 



(M-|-77P)+ • (M + ?/P) 



?7Q • K ■ Q + O (772), where the Hermitian matrix 
(M+ • M) ^, and the symmetrical tridiagonal ma- 



trix K has the elements 
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By using the correlation function (jH)) in the current 
density expression we can eventually write the fluc- 
tuation Hall conductivity in the form: 



3 



Aaxy 



hs f 



e- n I - duj' r dq' ^ , ^ 



n=0 



AT, 



ri'=0 



+ 2/ V^' + 1 £„'+ 1 (Qri+l,ri'Qri' + l,n + Q)i+l,?7' + lQn',Ti) 

I 



n'=0 



(9) 



where we have exphcitly specified the cut-off in the Lan- 
dau level sum. We note also that the zeroth order in r\ 
gives no contribution to Aoxy , as also established on the 
general grounds of the Onsager relations!^ The electric 
field enters Eq. © through the parameter /, defined in 
Eqs. (O- In order to apply the expression for Aa^y also 
in the limit i? — > 0, one has to expand the Q-matrix el- 
ements up to the linear term in /, namely Q = Q'-'^-' -I- 
/Q(i) O (p). Since Q(o) is diagonal and has the ele- 
ments Q^m — {w'^ + e^i^/} ^ ; one needs from Q'^' only 



the elements Qi+i_„ = 2y/n + lQrmQ^7i,n+i ^ 
and finally obtains in the linear response approximation; 
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which matches the formula found in Ref. 8 and is also 
coincident, after performing the g'-integral and replac- 
ing 77 = ~{kBT/eF)ot, with the expression given by 
Rcf. 9. If one neglects the cut-off procedure (i.e. makes 
Nc — + 00), and also the renormalization (i.e. replaces e 
by e), the result pO|) agrees with the Hall conductivity 
expression in weak electric field and arbitrary ma gnet ic 
field found from microscopical calculations in Ref. l2(]|i^ 
Taking further also the limit of weak magnetic field, one 
can find the expression given by Ref. |^ in the framework 
of the Boltzmann kinetic equation. 

Starting from Eq. JHJ we are able to compute also 
the fiuctuation Cooper pair density ^jV'l^^, keeping only 
the dominant zero-th order term in 77, and write the self- 
consistent equation for the renormalized reduced temper- 
ature parameter e, as already found in Ref. 0: 



e = In ■ 



T 



(11) 

Here, we have introduced the parameter g = 
(ttS-^s) and we have taken into account 
the expression of the quartic term coefficient^ b = 
2^oK'^e'^h? /m^ , with k the Ginzburg-Landau parameter 
K = Aq/^o- In analogy with the Gaussian fluctuation 
case, we shall adopt as definition for the critical temper- 
ature Tc{E, B) the vanishing of the reduced temperature, 
e = 0. The relationship between Tq and Tc(0, 0) = T^o, 
corresponding to Eq. Hll|) taken in the zero-fields limit 
at e = 0, has been already found in Ref. 5 and writes 
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We shall take as example the optimally doped 
YBa2Cu307_2:, for which typical characteristic param- 
eters are: s = 1.17 nm, ^0 = 1-2 nm, ^oc = 0.14 nm, 
K = 70 and Tco = 92 K. We assume a positive (hole-like) 
normal-state Hall conductivity a^y obeying the Ander- 



son's formula cF^^I'^xy = AT'^ with a generic A = 0.07 
at B = 1 T, and a linear extrapolation for the nor- 
mal state resistivity vanishing at T = 0, with a typical 
value pf^ = 84 ^f]cm at T = 200 K. The cut-off parame- 
ter c = 1. For the Fermi energy we take = fcsTO^Kjii 
while the parameter a will be given the positive value 
a — O.OlfiSiS^ in order to have 77 < and Aaxy < 0, and 
thus to account for the Hall effect's sign change occuring 
in the transition region. 

It should be mentioned that the relevance of the Af'- 
sign to the sign change of the Hall effect is still open 
to debate. The conventional s-wave weak coupling BCS 
theory predicts a positive excess Hall effect in the under- 
doped cuprates for a hole-like Fermi surface, in contrast 
to the experimental reportSf^^ although recent theoretical 
approaches based on the presence of preformed pairs 
on an additional contribution to the particle-hole asym- 
metry coming from the quadratic electron spectrumS^ or 
on the proximity of an electronic topological transitionSSi 
point out the possibility of an electron-like Hall sign in 
the hole-like doping range. Our purpose is however to 
illustrate the high electric field effect on the fluctuation 
Hall conductivity, for which a (or, equivalently, 77) rep- 
resents merely a prefactor that could be inferred from 
fits to the measurements. The experimental study of 
the non-Ohmic Hall effect could provide thus a supple- 
mentary and better tool for assessing the particle-hole 
asymmetry parameter a, without the uncertainty intro- 
duced by the previousiS*2S need to estimate the back- 
ground normal state contribution a^y, because the dif- 
ference a„j (E) - (Jxy (0) = Aaxy (E) - Aaxy (0) would 
be independent of a^y. An eventual experimental evi- 
dence of the non-Ohmic Hall effect behavior would also 
bring a strong argument in favour of the superconducting 
fiuctuations in the long-lasting debate on the causes of 
the Hall effect sign change in HTSC. 

Figure ^ shows the Hall conductivity axy — axy + 
Aaxy, normalized to the magnetic field value, for i3 = 1 
T (dotted lines) and B = 2 T (solid lines) at different 
magnitudes of the in-plane electric field. One can no- 
tice a strong suppression of the fiuctuation contribution 
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FIG. 1; Hall conductivity normalized to the magnetic field _B as a function of temperature for two values of the magnetic field, 
at several magnitudes of the (a) in-plane electric field; or (b) current density. A detail is shown in the inset. 



to the Hall conductivity ^a^y due to the high electric 
field. A higher magnetic field leads to a stronger reduc- 
tion of the fiuctuation Hall conductivity per se and thus, 
the fiuctuation suppression as a function of the electric 
field becomes relatively smaller. For E > 400 V/cm the 
fiuctuation part ^a^y becomes negligible with respect to 
the normal-state component cr^, so that the Hall con- 
ductivity turns out to be much more sensitive to high 
electric fields than the longitudinal one a^x- The latter 
preserves a significant superconducting fiuctuation com- 
ponent even for fields E = 1500 V/cni and B = 11 T, as 
shown in Ref. 0. This might be not too surprising, since 
Oxy is altered much more than Oxx also by a magnetic 
field alone (Ref. El). 

Throughout an experiment, constant current density j 
can be achieved much easier than constant E. Thus, 
Fig. QJd presents the same non-Ohmic effect on Oxy 
with j as the parameter. For this purpose, the equa- 
tion 2 = {oxx (T) + Aaxx (T, E, B)) E was firstly solved 
with respect to E at fixed T, B and j, where for the 
longitudinal fiuctuation conductivity A.axx the result of 
Ref. ^ has been used: 

e2 h r d^' r dq' . 

fis / J-oo 27r 27r ^ 

(12) 

One can see in Fig. ^ that, due to the very low re- 
sistivity at lower temperatures, high electric field values 
are difficult to attain, so that the non-Ohmic effect can 
be discerned only at the beginning of the transition, and 



only for current densities higher than 2 MAcm^^. This 
explains why the non-Ohmic effect on the Hall conduc- 
tivity has not been detected experimentally so far. Nev- 
ertheless, attaining current densities of a few MAcm"^, 
with minimal self-heating, on cuprate thin films of a typ- 
ical d = 100 nm thickness, might not be such a diffi- 
cult task, if one uses very short current pulses of the 
order of tens of nanoseconds, so that only the phonon 
mismatch at the film-substrate interface practically con- 
tributes to the sample temperature rise>^ According to 
literature data,^®'^^ a thermal boundary resistance of 
about Rhd — 0.5 mK-cm^/W between YBa2Cu307_2; and 
the substrate MgO would imply at about 100 K (where 
~ 40 pficm) a temperature rise AT — dRbdP Pxx — 
0.8 K for a current density j — 2 MAcm~^. This accuracy 
could be sufficient to discern the non-Ohmic behavior. 

In summary, we have calculated the critical fluctua- 
tion Hall conductivity for a layered superconductor in an 
arbitrary in-plane electric field and perpendicular mag- 
netic field in the frame of the TDGL theory using the 
self-consistent Hartree approximation. The main result 
is the formula that was found to reduce to previous 
results in the linear response limit Eq. (|10|l . Qualita- 
tively, high electric fields result in a strong suppression 
of the fiuctuation contribution to the Hall conductivity, 
in particular in moderate magnetic fields where order- 
parameter fluctuations are still strong. 
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